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We calculate the total cross section for single charged slepton production in association with a 
top quark at hadron colliders in the baryon triality (B3) supersymmetric model. We compute event 
rates for the Tevatron and LHC. We study the signatures for different supersymmetric scenarios 
including neutralino and stau LSPs. We perform a detailed analysis with basic cuts for the B3 
operator A231 using Monte Carlo simulations to show that the signal can be distinguished from the 
background at the LHC. In particular we employ the resulting lepton charge asymmetry. 



I. INTRODUCTION 

Supersymmetry (SUSY) [l| is a widely considered ex- 
tension of the Standard Model (SM) of particle physics 
0. If it exists, it is necessarily broken, with a mass 
scale of order of the TcV energy scale This en- 
ergy region is probed at both the Tevatron and LHC; 
the search for SUSY is therefore of paramount inter- 
est [1, [1] . It is the purpose of this paper to consider 
a specific supersymmetric production mechanism and 
investigate its viability at the Tevatron and LHC. 

The general renormalisable superpotential with mini- 
mal particle content [1] includes the following lepton 
or baryon number violating interactions. 
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-eat^'LtH', + U,yA'l,kU-D]Dl , (I.l) 



where we have employed the standard notation of 
Ref. for the superfields, couplings and indices. If all 
terms are simultaneously present, they lead to rapid 
proton decay Q . Therefore SUSY, must be augmented 
by an additional symmetry. The discrete anomaly-free 
gauge symmetries R-parity @ and proton hexality, Pq 
\W^ . forbid all of the above terms. However, R-parity 
does not forbid the dangerous dimension-five proton 
decay operators 

An equally well motivated solution to the proton decay 
problem is baryon triality, B3, a discrete anomaly- free 
^^-symmetry, which prohibits the UDD operator in 
Eq. dLH) [m. El, m . This solution has an additional 
feature, it naturally leads to small neutrino masses 
[13, [S [il, liZi, il^, ilB, i20|] , as experimentally observed 
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[2lj. Furthermore, B3 supersymmetric models which 
include also a candidate for dark matter can for exam- 
ple be reahzed within the UMSSM [2^ . 

The baryon-triality collider phenomenology has three 
main distinguishing features, compared to the Pg case 
0,111: 

1. the lightest supersymmetric particle (LSP) is not 
stable and can decay in the detector. It also need 
not be the lightest neutralino; 

2. SUSY particles are also produced singly, possibly 
on resonance; 

3. lepton flavour and number are violated. 

These lead to dramatically different signatures at 
hadron colliders [2^, [H, H^, HI] compared to the more 
widely studied Pg case. 

In the following, we focus solely on the signatures due 
to a non-vanishing LiQjDj- operator. At hadron col- 
liders this allows resonant single slepton and sneutrino 
production via incoming quarks, 



Uj + dk > I, , 

dj + dk -—^ Vi ■ 



{1.2) 
(L3) 



Here: Uj and dk denote up and down type quarks of 
generations j and fc, respectively; a bar denotes an 
anti-quark; £~ and Vi denote negatively charged slcp- 
tons and sneutrinos of generation i, respectively. 



The tree- level processes Eqs. p.2|) . (jl.3[) were first con- 
sidered in Refs. [13, HI]. Like-sign dilepton events or 
three lepton final states were considered in Refs. [1^, 
[sol , [sij . These papers assume a neutralino LSP, which 
can decay leptonically via j,, e.g. 
The case of a gravitino LSP was considered in Ref. [32] . 
The process has also been studied by the DO collabora- 
tion at the Tevatron [Hjlsj, for the operator L2Q1D1 
and a neutralino LSP, setting limits on the relevant 
masses and couplings. In Rcfs. [H, [sl, [13, [11] the 
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0.019 X (mj^/lOOGeV) 
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0.28 X (mj^/lOOGeV) 


•^231 
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-^232 
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-^331 


0.45 (m, = 100 GeV) 


-^332 


0.45 (mg = 100 GeV) 


•^i33 


0(10-") 







TABLE I: Upper 2a bounds on X'isk- The strong bounds on 
Ai33 stem from neutrino masses rriv, assuming < leV 
and left right mixing in the sbottom sector. The limits 
depend on the squark masses, rriq^^^^ is the mass of the 
left (right) handed squark qL(R)- 



£^ , respectively. The main production mechanisms for 
production, at O(a^), are 

g + g t + t J 

where q (q) is a (anti-)quark. This is only kinemati- 
cally allowed if 




FIG. 1: Feynman diagrams contributing to the partonic 
process g + dk ^ t + £~ . 

NLO-QCD corrections were considered and in Rcf. 
jsrl l the SUSY-QCD corrections were taken into ac- 
count. Gluon fusion contributions were included in 
Ref. [13. 

The case j = 3 in Eq. (jl.2p is special, as there are 
no top quarks in the incoming proton. Instead, one 
must consider the production of a single slepton in 
association with a SM particle. Several mechanisms 
for associated single supersymmetry production, e.g. 
djdk Xi^T ^ have been studied in the literature, see 
for example Refs. [Ill, M, S SI, (H] . In the follow- 
ing, we investigate in detail the case of the operator 
LiQ^Dk- Here, single charged slepton production is 
only possible in association with a top quark. Be- 
fore studying the phenomenological details, we first 
recall the strongest experimental bounds on the cou- 
plings A;c,t. at the 2a level. They are shown in Table [J 
0, m, Ea S, [13 ■ We neglect bounds, which assume a 
specific (Standard Model) quark mixing between the 
three generations (48j or bounds using the rcnormal- 
ization group running of A -gj, [H, [13, \^ . 

At leading order there arc two production mechanisms 
for slepton production in association with a top quark. 
First the Compton-like processes 

g + dk i^+t, (I.4a) 
ff + rffe -> it+i. (I.4b) 

The relevant leading-order diagrams are given in 
Fig. [T] Here, g denotes an incoming gluon in the pro- 
ton and t a final-state top quark. 

The second slepton production mechanism is tt pair 
production followed by the t or t decaying into if or 



mt > TOfj -t- nid^ . (1.6) 



Since, as we shall see, the branching fraction for the 
B3 top quark decay is small, we only consider one B3 
decay, for either the top or the anti-top quark. 

In Ref. [3§| , single slepton production was considered 
for the specific case of A333 ^ 0. This process is how- 
ever disfavoured due to the strict bound on the rele- 
vant coupling from neutrino masses, cf. Table |T] [ssj . 
Thus the work was extended to the couplings X'^^i and 
^^332 ■ We go beyond this work to include a signal 
over background analysis. We also present the analytic 
formula for the cross section, Eq. (|II.9[) . for the first 
time, and analyse the resulting signatures. We give a 
detailed phenomenological analysis for the special case 
-^231 which can be generalized to X'i;ji. 

In Ref. [sot i top quark pair production and subsequent 
top decay via A^gj, was considered. Off-shell top quark 
effects were also taken into account. A signal over 
background analysis was performed for two scenarios. 
The first scenario assumed maximal stop-scharm mix- 
ing. It was pointed out that associated slepton produc- 
tion with slepton masses 150 GeV and 200 GeV can be 
measured, depending on the magnitude of A^g^,. The 
second scenario assumed no flavour violation in the 
squark sector. Ref. [5^ claimed that in this regime 
sleptons with mass 200 GeV can not be measured. 
We go beyond the work of [5^. We show that it is 
possible to detect associated slepton production even 
for slepton masses larger than 300 GeV, if or X[r^i 
is of 0(0.1). We will achieve this with the help of the 
Compton-like process p.4p . 

The outline of this paper is as follows. In Sect. [H] 
we calculate the cross section for the production of 
a charged slepton in association with a top quark, 
at leading order. In Sect. IIIII we systematically 
present the possible resulting signatures at the LHC. 
In Sect. IIVI we discuss in detail a case study for the 
operator X231L2Q3D1. We study the dominant tt and 
backgrounds. Using the HERWIG Monte Carlo 
program [5l|, [H, , we devise a set of cuts in order 
to distinguish the two. We do not include a simula- 
tion of the detector. Our conclusions are presented in 
Sect. El 
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II. SINGLE SLEPTON PRODUCTION VIA 

A. Partonic cross sections 

The spin and colour averaged matrix element squared 
for the Compton-like process Eq. (jl.4p is given at lead- 
ing order by 

(m? - m2)(m? - s - t) - (Bm^ - m? + s){i - m? ) 



{t - ml) 



2\2 



2K-s+ft-m|)(m| 
s{i-ml) 



(11.7) 



where as is the QCD coupling constant, Cp = 4/3 is 
the quadratic Casimir of SU{Z)c, rrij;^ is the mass of 
the slcpton and i^J^ is the relevant matrix element of 
the left-right slcpton mixing matrix. The explicit form 
as a function of the mixing angle is given, for example, 
in Rcf. [53|. In accordance with the parton model, we 
have neglected the mass of dk- We have made use of 
the partonic Mandelstam variables 

s = (dk + gf = (t + £^)^ , (II.8a) 
i^{dk-i^f^{g-tf, (II.8b) 

where we denote the particle four momenta by the 
particle letter. Integrating over phase space, we obtain 
the total partonic cross section: 



64s2 I 2s 



+ (J + 2m^)ln ^ + 



p\ 2ml{ml-ml){t+~t^) 



+ 



P+J P+P- 



2(mf^. +n4~ 2mlm]. - mis) ^ / 



2(t+-t_)(m? -to2_5) 



where 



P±=m^- t±, 



(11.9) 



(11.10) 



i± = m\ - 2 + "^1 - T (s, m|, to?)], 

(11.11) 

with the phase-space function given by \{x,y,z) = 
a;^ + 2/^ + 2^ — 2xy — 2xz - 2yz. 

The tree-level partonic matrix element squared for top 
quark pair production is given for example in Ref . [54| . 
We shall only consider on-shell top quark pair produc- 
tion. The slepton then arises through the decay of a 



real top quark. In order to obtain the signal rate, we 
thus also require the partial decay width of the top 
quark, via the LiQ^Dk operator. It is given by 



327rTOt 



2 2 2 \ 



(11.12) 



See also Refs. [39|, |48|, [SJ, |55|, Isg . We obtain a branch- 
ing ratio of 8.2 x 10~^ for the R-Parity violating top 
decay pLT2|) for X'-^^. = 0.1, mt = 175 GcV, top width 
= 1.5 GeV and m^^ = 150 GeV. We neglect the 
mass of dk and set L^'^ = 1. 



Total Hadronic cross section 



10^ 



10 



10 



10 



10" 



10'' 



E 1 ' 


1 ' 1 ' 1 


' 1 ' E 


- V 


pp^ 


£~t via A -31 : 




--■ PP^ 


ijt via A -32 - 


\ 


PP^ 


irt via A;33^ 


\ 
ft 

A 
i\ 

— ii 
= ) 
) 
I 
( 




= 1.96 TeV 


_ 


\ \ 
\ \ 
\ \ 
\ \ 

N \ 






\ \ 
*. \ ^ 
\ 

\ 




1 . 


\ 

\ 

\ 

N 

\ 

\ 

\ 

1 , 'i\ , N 


, 1 ,\ 



100 200 300 400 500 600 

m« [GeV] 



FIG. 2: Single slepton production in association with a top 
at the Tevatron. The cross sections for ift production are 
equal to the cross sections for £~t production. 

In Fig. [2] (Fig. [3]) , we show the hadron level cross sec- 
tion at the Tevatron (LHC) for single slepton produc- 
tion in association with a top quark, as a function of 
the slepton mass including both production mecha- 
nisms. We set A^3j, =0.1 and assume it is the only 
non-vanishing B3 coupling. We vary the index k and 
the charge of the final state slepton, which correspond 
to different parton density functions (PDFs) . Here we 
use the CTEQ6L1 PDFs [13 . The renormalisation, 
fiR, and factorisation, /ip, scales are taken to be equal, 
PR = Pf = rn, where to = 2TOt [= toj^ -|- mt] in the 
case of slepton production via a tt pair p.5p [via the 
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FIG. 3: Same as Fig.O but for the LHC. The cross section 
for ift production via Ai32 {K33) is equal to the cross sec- 
tion for £^t production via Ai32 {X'iii), as it always involves 
incoming sea quarks. 



Compton-likc process (jl.4[) ]. Furthermore, we set the 
left-right slepton mixing matrix element i^'^ equal to 
one. Results for other values of A^g^. and mixing matrix 
elements L^'^ are easily obtained by rescaling accord- 
ing to Eqs. pi.gp and pi.l2[) . The top mass is taken 
to be 175 GcV and the total (SM) top quark decay 
width to be 1.5 GeV. We take m^g = rrib = 4.5 GeV, 
if we have a b quark in the final state and neglect the 
masses of the d and s quarks. 

In both figures, we see a kink in the cross section 
when m^^ = rrit — rriiif, . For smaller slepton masses the 
top quark pair production mechanism dominates; for 
larger masses the Compton-like processes dominate, 
since the slepton can no longer be produced on-shell 
in top decay. 

For comparative discussions later. Fig. [4] (Fig. [5]) shows 
the NLO hadronic cross section for resonant sneutrino 
production, cf. Eq. p.3[) . at the Tevatron (LHC) via 
A^3fe =0.1, including NLO QCD corrections [s^. We 
employ the MS renormalisation scheme and the (NLO) 
CTEQ6M PDFs The renormalisation and fac- 

torisation scales are taken to be the sneutrino mass. 

In Fig. [2l we see that at the Tevatron, even for small 
slepton masses, m^^ = 100 GeV, we expect only 25 (25) 
charged slepton events with negative (positive) charge, 
i.e. i~ (if), for an integrated luminosity of 1 fb~^ and 
the (relatively large) couphng A-g^ =0.1. The cross 



FIG. 4: Single sneutrino production cross section at the 
Tevatron. The cross sections for u* production are equal 
to the cross sections for Ui production. 
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FIG. 5: Same as Fig. HJ but for the LHC. The cross section 
for V* production via Ai32 (Ai33) is equal to the cross section 
for Vi production via Ai32 (Ai33), since only initial-state sea 
quarks are involved. 



section is dominated by the tt pair production (jl.5p . 
Only 10% of the above sleptons at the Tevatron are 
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FIG. 6: Factorisation scale = f-m and renormalisation 
scale IJ.R = r-m dependence of the hadronic £~t production 
cross section via A^ai at the Tevatron. /if and fiR are 
independently taken equal to 2 and 0.5 times m, where 
m = 2mt (= m^-+mt) in the case of slepton production via 
a tt pair Eq. (|I.5|I (via the Compton-like process Eq. ()I.4|) ) . 
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FIG. 7: Same as Fig. [5] but for the LHC. 



produced by the Compton-like process (|L4|. At the 
Tevatron. the cross section is symmetric in the slepton 
charge due to the charge symmetry of the incoming 
state. 

As we can see in Fig. [31 we have a significantly larger 
hadronic cross section at the LHC for a given slep- 
ton mass. In particular, for m^- = 100 GeV and 
A^3^ = 0.1 the LHC wiU produce more than 31 000 (26 
000) sleptons £^ {£^) for an integrated luminosity of 
10fb~\ 27% (11%) of these sleptons are produced via 
the Compton-like process. For the same coupling and 
for nij}. = 100 GeV, we will produce approximately 
14 000 sneutrinos at the Tevatron (Fig. [J) for lfb~^ 

and 3 800 000 at the LHC (Fig.O for 10 fb^^, via the which they show, for nif < rrit — nib, to be equiva- 



partonic process Eq. p.3p . Thus, depending on the de- 
cays, we might expect this to be the discovery mode, 
for equal supersymmetric masses. Here we focus on 
the potential of the charged slepton production cross 
section. 

For heavier charged sleptons, m^- — 800 GeV, we ex- 
pect no events at the Tevatron and more than 110 
(25) i~ (if) events at the LHC with 10fb"\ Above 
the threshold of m^^ = rrit — , practically all slep- 
ton events are produced via the Compton-like pro- 
cess, since the other process only proceeds via off- 
shell top quarks. The cross section is so small be- 
cause the parton luminosity is too small at the required 
high values of the proton/anti-proton fractional mo- 
menta, X > 0.1. This situation changes at the LHC, 
where we probe significantly smaller values, x < 0.1, 
for the same slepton mass. Furthermore, the Teva- 
tron will produce no sneutrinos, for X'^^^ = 0.1, and 
mpi = 800 GeV. For the same set of B3 parame- 
ters, the LHC will produce about 3 200 sneutrinos for 

At the LHC, there is an asymmetry between the 
hadronic cross sections for and £^ production via 
the LiQ^Di operator (fc = l!). This is perhaps not 
surprising, as the initial state is asymmetric under 
charge reversal. In the case of the Compton-like pro- 
cess (|I.4p . the asymmetry is due to the negatively 
charged slepton being produced by an incoming va- 
lence d-quark, while the positively charged slepton is 
produced by a d sea quark. The latter has a lower 
luminosity in the proton. In Sect. IIVI we will use this 
asymmetry to separate the B3 process from the SM 
background. 

In order to estimate the influence of higher order 
corrections on the production cross section, we vary 
the renormalisation and factorisation scales indepen- 
dently between m/2 and 2m. At the Tevatron, Fig. [SI 
the hadronic cross section for i~t production via X'^^^ 
changes by up to 40%. At the LHC, Fig. [3 the scale 
uncertainties are reduced to approximately 25%. In 
the domain where m^^ < rrit — nid^ , we have a stronger 
dependence on the renormalisation scale compared to 
> rrit — rrid^, because tt production is 0{al{^r))- 
According to Ref. [s^, [g^l, NLO-QCD corrections, in- 
cluding a NLL resummation, increases the tt produc- 
tion cross section by approximately 40% (80%) at the 
Tevatron (LHC). 

The hadronic cross section for single stau, f , produc- 
tion via a non- vanishing A333 coupling, was also con- 
sidered by Borzumati et. al. [3^. There, the 2 — > 2 
processes, Eqs. (|L4p were included, together with the 
(tree-level) 2 — > 3 slepton-strahlung processes 




t + b + r 



(11.13) 
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lent to the 2 — > 2 processes, Eqs. p.5l) . The b and 
f ^ are produced via a virtual top. They employed the 
CTEQ4L [6l| PDFs and all matrix elements were mul- 
tiplied by the CKM factor Vtb- We have calculated the 
hadronic cross sections using the same PDFs and the 
same parameter set [62j . We coincide exactly, where 
single slepton production is dominated by the tt pro- 
cess, i.e for TTif < rat — TUb- For rrif > rrit — mt, we 
underestimate the total cross section at the Tevatron 
by 20% for ruf = 300 GeV and by a factor of roughly 
two for iTLf = 200 GcV; compared to Ref. [3^. In this 
region the above 2-^3 processes, where the slepton is 
produced by a quark-antiquark pair, can give the main 
contribution compared to the gb — > ft partonic pro- 
cess, where a gluon and/or sea-quark is needed with 
large Bjorken x. However, in this region where there 
are large discrepancies practically no sleptons are pro- 
duced at the Tevatron. Our prediction for the LHC 
differs by -1-30% for mf > nit — fnt- 

Borzumati et. al. extended their analysis to the 
A332 and A331 couphngs [4l[. They present the re- 
sults for the 2-^2 process (|I.4p and the 2-^3 
process pi.l3[) . separately. For nrif < rrit — rnd^, we 
coincide exactly at the Tevatron as well as at the LHC. 
For rrif > rrit — ^dj. , our predictions coincide exactly 
with their cross section predictions for the 2^2 
process. Furthermore, it is shown in [4l[ that for 
TOf > rrit — n^dk , the 2 — s- 3 contributions are small 
or even negligible. At the Tevatron, the 2^3 process 
contributes roughly 35% (5%) to the total hadronic 
cross section for A332 ^ (A332 7^ 0). At the LHC 
these contributions are 25% (5%). The reason is that 
the cross sections induced by the 2 — > 3 process have 
similar sizes for any value of k. But the 2—^2 pro- 
cess for A332 7^ (A33]^ 0) is enhanced by a factor 
of 5 ( > 10) due to a s-quark (valence d-quark) in the 
initial state. 

We conclude, that our LO approximation is valid in 
the phenomenologically relevant region, where one is 
able to produce a single slepton in association with a 
top quark. We have not included the 2 — » 3 processes 
as they are formally higher order. Furthermore, the 
essential ingredient in our phenomenological analysis 
below is the lepton charge asymmetry due to a non- 
vanishing A^3]^ coupling. The 2^3 processes do not 
contribute, as their initial states are charge symmetric 
and their contributions to the hadronic cross section 
are only 5%. 

We end this section by presenting in Table |TT] selected 
cross section predictions for slepton production with 
= 100 GeV, m^^ = 250 GeV and m^^ = 800 GeV 
at the Tevatron and the LHC via X'^^^ =0.1. 





Tevatron 


LHC 


m|- = 100 GeV 


25.5 fb 


3180 (2620) fb 


m^-- = 250 GeV 


2.10 X 10"^ fb 


259 (80.0) fb 


m^-- = 800 GeV 


2.86 X 10"^ fb 


11.6 (2.54) fb 



TABLE IL Hadronic Cross section predictions for t (ift) 
production via Ksi = 0.1 at the Tevatron (TS = 1.96 
GeV) and the LHC 0/5 = 14 GeV). Results are presented 
for the CTEQ6L1 [sl PDF parametrisation. 



ing mixing between left- and right-handed sleptons the 
possible tree-level decays are: 




' bdk 

' Xn 



(IIL14) 



The branching ratios depend on the masses of the spar- 
ticles, the admixtures of the gauginos and on the size 
of the A^3j. couphng. We shall first assume, that the 
lightest neutralino, Xii is the lightest supersymmet- 
ric particle (LSP). Possible decay modes via the A -3^ 
interaction are: 



-0 >■' 
Xi — ' 



Ctdk 




(IIL15) 



Here the branching ratios depend mainly on the ad- 
mixture of the lightest neutralino. The heavier neu- 
tralinos X2 3 4 ^-^^^ the charginos X12 dominantly decay 
into lighter gauginos via gauge interactions, as in the 
Pg-MSSM. We have neglected the decay Xi ^ 
which is suppressed except for very light neutralino 
masses 



In SUSY scenarios, where the slepton (sneutrino) mass 
is of the order of a few hundred GeV, the slepton 
(sneutrino) will decay dominantly into the lightest 
neutralino and a lepton (neutrino). However, sig- 
nificant chargino decay modes are also possible, if 
they are kinematically accessible. Furthermore, de- 
cay chains involving a top quark in the final state are 
either phase-space suppressed or kinematically forbid- 
den, unless the slepton is very heavy. This affects 
the neutralino decays piI.15[) involving charged lep- 
tons. Therefore, the dominant hadron collider signa- 
tures of single slepton production in association with 
a top quark are 



(IH.16) 




III. POSSIBLE LHC SIGNATURES 

Apart from the B3 process, the sleptons and sneu- 
trinos can decay through gauge interactions. Neglect- 



In parentheses are the neutralino LSP decay products 
piL15P : the particles in brackets arise from the top 
quark decay. As mentioned before, for fc = 1 there is 
an asymmetry between the number of positively and 
negatively charged leptons if at the LHC. 
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The dominant signatures for a resonantly produced 
single sneutrino are 

[ bdk 

hdk^h ^ <^ u,{D,hdu) , (111.17) 
[ i^i {vibdk) 

again the neutralino decay products are in parenthe- 
ses. Although the sneutrino production cross section 
at the LHC (Fig. [5]) is up to two orders of magnitude 
larger than the slepton plus top quark cross section 
(Fig. [3]), the event signature (jIII.17[) is much harder to 
extract above the SM background. It involves only two 
jets and possibly some missing transverse energy. It 
therefore suffers from a large QCD background. How- 
ever, if the sneutrino decays into charginos and heavier 
neutralinos are possible (jIII.14p . we can have (addi- 
tional) charged leptons in the final state. 

We now consider SUSY scenarios, where the scalar tau 
(stau) is the LSP instead of the lightest neutralino 
0, HI, [fi^. In this scenario the lightest neutralino 
dominantly decays into a tau and the stau LSP, — > 
fj^r^. For i = 1,2, the stau will dominantly decay 
into a tau and a virtual neutralino, leading to a four- 
body decay of the stau LSP [1^ . The signatures for a 
stau LSP are 

r €-T± ij^vibdu) [hw+] 

gdk^^-t^l . (in.18) 

i ^-T± {T^vMk) [hW+] 

The particles in parentheses are now the stau LSP 
decay products and the particles in brackets are from 
the top quark decay. The difference between the final 
states in Eqs. ()III.18jl and (jlll.iep is, that for a stau 
LSP, the event is accompanied by an additional pair 
of taus compared to scenarios with a neutralino LSP. 
We find the same behaviour for the sneutrino decay 
chains. It is therefore easier to distinguish the signal 
from the background in stau LSP scenarios as long 
as one is able to reconstruct the tau pair in the final 
state. 

Note that for i = 3 the two-body stau decay is kine- 
matically suppressed, or forbidden, due to the large 
top quark mass. The stau LSP will in this case decay 
via a virtual top quark. Furthermore, we can pro- 
duce heavy staus, t2, as well as light staus, fi, due to 
left-right mixing in the stau sector. In this case the 
signatures are 

{T-T+ (bdkW-) [hW+] 
T-T- {hdkW+)[bW+] , (III.19) 
Z^/h° {bduW-) [hW+] 

and 

gdk ^ fft ^ {hdkW-) [bW+] . (III.20) 



'^231 





0.1 


0.2 


0.3 


0.4 


Br(/i- ^t + d) 
Br(^- -^^^-+ xi) 
Bi ifl- ^ ^i- + X'i) 
Br(/i" i^^c + xl) 


0.0% 
90.9% 
3.2% 
5.9% 


2.2% 
88.9% 
3.1% 
5.8% 


8.4% 
83.3% 
2.9% 
5.4% 


17.1% 
75.4% 
2.6% 
4.9% 


26.8% 
66.5% 
2.3% 
4.3% 



TABLE III: Relevant branching ratios for SPSla' for dif- 
ferent couplings A231. 



The particles in parentheses are the stau LSP decay 
products and those in brackets are from the top quark 
decay. We see in Eq. piI.lQp that one of the f2 decay 
chains involves like-sign tau events. This can help to 
distinguish signal from background although poor tau 
identification could limit this possibility. 

IV. NUMERICAL STUDY FOR A231 / AND 
A xi-LSP 

A. The Scenario and Basic Cuts 

We now perform an explicit numerical study of sin- 
gle associated slepton production. We focus on the 
more difficult case of a neutralino LSP and restrict 
ourselves to A231 ^ 0, as the dominant B3 coupling. 
We assume that similar results can be obtained for 
X'l^i 7^ 0. A central analysis criterion will be the lep- 
ton charge asymmetry of the final state. 

According to Eq. piI.16|) . the final-state signature to 
examine is 

q + t + {b + d + iy) + [b + w'^] , (1V.21) 

with the decaying hadronically. We thus have 
one charged Icpton, some missing px, and five jets, 
where two are b-quark jets. In our specific scenario, 
the charged lepton is a muon. 

The main background for this process is tt + j pro- 
duction (which has recently been calculated at NLO 
[1^) followed by the semi-leptonic decay of one of the 
top quarks. The second background we examine is 
bb + + jets production followed by the leptonic 
decay of the W boson. 

For our simulation, we assume an SPSla' similar sce- 
nario [63]. We take the SPSla' spectrum and couplings 
and add one B3 coupling, A'231. The relevant SPSla' 
masses are: 

mi± = 190 GeV; m£>^=173GeV; (IV.22a) 
m^o = 98GcV; m^o = 184GeV; (IV.22b) 
m-± = 183 GeV. (IV.22c) 

All the charged slepton decays of Eq. (|III.14|) are 
therefore kinematically possible. The corresponding 
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branching ratios are given in Table [Till for various cou- 
plings A231 . Note that kincmatically the sneutrino can 
only decay via the neutralino or via the coupling. 
The potential signature would then be two jets possi- 
bly with some missing energy, cf. Eq. piI.17p . 

For the simulation of the single slepton plus top quark 
signal we have written our own Monte Carlo program 
using the Les Houches accord [t^ and linked this to 
HERWIG6.5 [m, m, [nj. The averaging of the colour 
flow in the s- and i-channel single slepton produc- 
tion diagrams is implemented by the method devel- 
oped in Ref. [t^ . The supcrsymmetric particle spectra 
are produced with SOFTSUSY The tt-background 
is simulated using the MCQNLO program [tJ, Its} . 
The bb + -I- jets background is simulated by using 
MadEvent [t^ to generate a sample of 66-1- + 2 jet 
events which are then showered and hadronizcd us- 
ing HERWIG6.5. We use the CTEQ61 parton distri- 
bution functions [s^l- The top quark mass is set to 
mt = 175 GeV. 



0.1 




50 100 150 200 250 
[GcV] 

FIG. 8: Relative pr distribution of the final-state signal £^ 
at the LHC for SPSla', employing only the isolation cut 
on the lepton. 



0.01 



0.001 



lc-04 




5(1 100 15(1 200 250 
[GeV] 

FIG. 9: Relative pr distribution of the final-state £^ from 
tt+j background at the LHC, employing only the isolation 
cut on the lepton. 

Since our signature is very similar to the final state 
and distributions of tt + j production followed by the 
semileptonic decay, we use the standard set of CMS 
cuts for tt production followed by the semileptonic de- 
cay, given in Ref. [t^] and require an additional jet. 



This set of cuts leaves the large semileptonic tt + j 
production, for which the cuts are designed, and fewer 
bb + + jets events as backgrounds for the signal 
process. The precise cuts are summarised below. 

The main difference between the semileptonically de- 
caying top pair and our signal is the pt distribution 
of the lepton stcming from the slepton compared to 
the one from the from one of the top decays. We 
therefore compare in Figs. [8] and Fig. [9] the px distri- 
butions of the leptons arising from the signal and the 
tt + j background processes, respectively. We see, that 
the Pt of the signal leptons has a peak around 50 GeV. 
This peak corresponds to the mass difference between 
the slepton and the neutralino with the energy car- 
ried away by the lepton subtracted. The background 
lepton distribution peaks at 25 GeV and then falls 
more steeply than the signal distribution for increas- 
ing Pt- We thus harden the CMS semileptonic tt cut 
for the isolated observed lepton from px > 20 GeV to 
PT > 35 GeV. 

In addition to the charged lepton in the final state, we 
require two tagged 6-jets, as well as two further jets. 
Thus, the employed cuts are 

• 1 isolated lepton with pseudo-rapidity 77 < 2.4, 
PT > 35 GeV. The isolation cut required less 
than 2 GeV of transverse energy in a cone of ra- 
dius 0.4 around the lepton direction. 

• 2 isolated 6 jets and 2 non-6 jets, pseudo-rapidity 
T] < 2.4, Pt > 30 GeV. 

The jets are defined using PXCDNE [6§| which uses the 
mid-point between two particles as a seed in addition 
to the particles themselves to improve the infrared be- 
haviour of the algorithm. A cone radius of 0.5 was used 
to define the jets. For the bottom and charm quarks 
produced in the perturbative stage of the event the 
nearest jet in {r],(j)) is considered to have been pro- 
duced by that quark if the distance in (r],(j)) was less 
than 0.2. We employ a 6-tagging probability of 0.6 
and the probability for mistagging a c-quark or light 
quark as a 6-quark of 0.05 and 0.02 respectively. 

For the signal, we simulated 10^ events. Employing 
all cuts, including pt(^*) > 35 GeV, we have 5 x 10'^ 
surviving £~ events and 1.7 x lO'^ surviving i'^ events. 
For the 66 -I- + jets background we simulated 10^ 
events for both and production. After all 
cuts we are left with 2.9 x 10'' £' and 3.0 x 10'' i+ 
events, respectively. 10^ tt + j events were simulated 
resulting in 1.35 x 10^ events for £~ production and 
1.36 X 10^ events for ^"''production. This is summarised 
in Table [IVl 

For the simulated signal, we set A231 = 0.053. In the 
following we will estimate the signal for other values 
of A231 by taking into account the A231 dependence of 
the cross section. We also employ the A231 dependence 
of the the p" — > /i~ -I- xl branching ratio. 
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simulated 


£ -events after cuts 


-events after cuts 


Events -pb ^ 


signal 


99 900 


5 042 


1664 


0.0108 


W^"+66+jets bg 


994 000 


28 600 





0.0431 


W+ + bb + ]ets bg 


993 500 





29 700 


0.0625 


tt + lj bg 


9 990 500 


135 330 


136 360 


22.00 



TABLE IV: Results of simulating SPSla' with cuts given in the text. The number of leptons and the expected event 
rates are after cuts. 



B. Lepton Charge Asymmetry 

In order to distinguish the signal from the background 
at the LHC after these cuts, we propose as the decisive 
observable the lepton charge asymmetry 



(IV.23) 



Here Ng+ and iVf_ are the number of events with a 
positively or negatively charged lepton, respectively. 
In Fig. [31 we can see the separate signal cross sec- 
tions for and £~ production at the LHC. For 
mf± > rrit — m^^ the £^ cross section is significantly 
larger. This is due to the fact that the c?-quark PDF 
luminosity is significantly larger than that of the d- 
quark for x > 10~^. 

We would expect the lepton charge asymmetry to be 
zero for the ti + j background, as we have an equal 
number of top quarks and anti-top quarks. For the 
background process bb + + jets, we expect a posi- 
tive asymmetry, since the (valence) u-quark luminosity 
is significantly larger than the (valence) d-quark lumi- 
nosity in the proton. For the signal, as we have seen, 
we expect a negative asymmetry. 

However, in general, inclusive tt production has a 
charge asymmetry in the final state at the LHC. It 
has been shown to be in the range [—0.025%; 0], if the 
detector has a symmetric acceptance in the rapidity 
range [—yo'-,yo]- For i/n — > oo (0) the asymmetry goes 
to (-0.025%) [zl Ei. This stems from the asymme- 
try in qq induced tt production, which in turn is due 
to the interference of C-odd and C-even modes, where 
C is the charge conjugation operator. In the follow- 
ing, we will neglect this small asymmetry because the 
statistical fluctuations lead to an even larger asym- 
metry. The number of £^ events in Table IIVI for the 
tt+i background are consistent with a lepton charge 
asymmetry of zero within two sigma. 

In Fig. [ini we show the significance, E, of the signal 
for the SPSla spectrum as a function of A231, where 



{A 



SM 



A 



SM+S) 



AA 



SM 



(IV.24) 



background combined. A Asm is the error of the SM 
asymmetry prediction assuming purely statistical er- 
rors for the number of positive and negative charged 
leptons for each process separately, i.e. y/Ni^ and 
Ne-. The significance is shown for integrated lumi- 
nosities at the LHC of 30 fb~^ (lower curves), 100 fb~^, 
300fb"\ and 1000 fb~\ respectively. We vary the 
cross section by ±20% (grey region) to show possible 
effects due to higher order corrections for the signal 

(c/. Fig. ED . 




FIG. 10: Significance at the LHC as a function of A231 
for SPSla' with lepton pr > 35 GeV. We show the signifi- 
cance for an integrated luminosity of 30 fb^^ (lower curve), 
100fb-\ 300fb-\ and 1000fb-\ respectively. Further- 
more, we varied the signal cross section by ±20% (grey 
region) . 



In Fig. [ini we see that for 30 fb^^ we can probe cou- 
plings down to about 0.3 for the SPSla' spectrum. 
In the SPSla' spectrum the squark mass is 544 GeV, 
thus the experimental bound is A231 < 1.0, cf. TablelH 
For 300 fb~^ we can probe couplings down to about 
0.15. In the extreme case of 1000 fb^^ this improves 



to about Ao 



231 



0.1. 



Here Asm is the SM lepton charge asymmetry. 
AsM+s is the asymmetry for the signal and the SM 



We have repeated the above analysis for the parameter 



-^231 





0.1 


0.2 


0.3 


0.4 


Br(/i' t + d) 
BTifi- ^ 11- +X\) 
Br(/i- -^^i- + X2) 
Br(/i~ +Xi") 


0.0% 
60.9% 
13.8% 
25.3% 


22.0% 
47.5% 
10.8% 
19.7% 


53.0% 
28.6% 
6.5% 
11.9% 


71.8% 
17.2% 
3.9% 
7.1% 


81.9% 
11.0% 
2.5% 
4.6% 



TABLE V: Relevant branching ratios for SPSlb for differ- 
ent couplings A231. 



set SPSlb [bI]. Here we have the following masses: 

m^^ ^ 342 GeV; m£>^, 333 GeV; (IV.25a) 
m^o = 163 GeV; m^o= 306GeV; (IV.25b) 
m-± = 306 GeV. (IV.25c) 

We show the branching ratios for different A23J in Ta- 
ble |Vl We see, that the B3 decay into a d quark and 
a top quark is the dominant decay for large A231, i.e. 
A231 > 0.19. One might thus consider an analysis 
based on this decay mode. However the signature is 
ti + j, which has a very large background. We thus 
continue to consider the neutralino decay mode. The 
significance will then approach a constant value for a 
constant luminosity and large A231, because the cross 
section and the B3 decay both scale with A231. Fur- 
thermore the slepton mass is now significantly larger, 
but so is the lightest neutralino mass. The mass differ- 
ence however has grown, leading to significantly higher 
charged leptonpr's compared to SP SI a', c/ Fig. [51 We 
thus impose the stricter cut on the lepton transverse 
momentum 

Pt(^=^) > 70GcV. (1V.26) 
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A23I 



FIG. 11; Same as for Fig. 1101 but for the parameter set 
SPSlb with lepton pr > 70 GeV. 

The results are shown in Fig. [TlJ In this case, for 
the relatively low integrated luminosity of 30 fb~^ we 
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-^231 





0.1 


0.2 


0.3 


0.4 


Br(/i" t + d) 
Biif,- -^^^-+ xi) 


0.0% 

100% 


23.7% 
76.3% 


55.4% 
44.6% 


73.6% 
26.4% 


83.2% 
16.8% 



TABLE VI: Relevant branching ratios for the high pr sce- 
nario for different couplings A231 . The scenario is described 
in the text. 



have no chance of observing the signal via the lepton 
asymmetry; the neutralino branching fraction is too 
small to have enough events. In fact, it is only for 
the extremely high integrated luminosity of 1000 fb~^ 
that we have a significant sensitivity range, down to 
about A231 = 0.2. 

In order to see what can be probed at the LHC, we 
have chosen as a third example a mass spectrum which 
optimizes our signal. For this we considered a modified 
SPSlb spectrum, where we first lowered the mass of 
the lightest neutralino to 

m^o = 80 GeV, (IV.27) 

in order to obtain a larger mass difference between 
the smuon and the lightest neutralino. We can then 
harden the pr cut to 

Pt(^=^) > 120GeV. (IV.28) 

This leads to a better signal to background ratio com- 
pared to SPSlb. Second, we increased the masses of 
X2 and xt to 

m^o=m~± = 450 GeV. (IV.29) 

This increases the jj^ M +Xi branching ratio com- 
pared to SPSlb, because decays into heavier neutrali- 
nos and into charginos are now kinematically forbid- 
den. We show the relevant branching ratios for differ- 
ent A231 in Table IVTl We refer to this scenario as the 
high-pr scenario. The resulting significance for the 
high-pT scenario is shown in Fig. [T^l 
As can be seen, for an integrated luminosity of 30 fb^^ 
we still have no sensitivity in A232 . But now for an inte- 
grated luminosity of 300 fb""'^, we can probe couplings 
down to 0.19, well below the experimental bound of 
1.5, cf. Table U where now n^^ ^ 830 GeV in SPSlb. 
For an integrated luminosity of 1000 fb~^ we can probe 
couplings down to 0.11. 

The influence of systematic errors in the background 
cross section on our sensitivity are small. Varying 
the tt -I- J cross section by -|-10% (—10%) changes 
the asymmetry by roughly —9% (-1-11%). Varying the 
bb + W"^ -|- jets cross section by ±10% only effects the 
asymmetry by ^1.6% for SPSla' with A231 = 0.3 and 
by =f1.2% for the high-p^ scenario with A231 = 0.3. 
Yet, detector effects resulting in an error on the ob- 
served charge asymmetry are a problem. Misalign- 
ment in the detector can lead to a difference in px 
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•^231 



FIG. 12: Same as for Fig. 1101 but for the high-pr scenario 
and with lepton pT > 120 GeV. The scenario is described 
in the text. 

measurement of positive and negative leptons, respec- 
tively. This will lead to an observed, effective charge 
asymmetry after a cut on the lepton pr [13] . An anal- 
ysis of this must be performed by the experimentalists 
and is well beyond the scope of this paper. 

For SPSla' with A23X = 0.3, a simulated detector 
based charge asymmetry of 0.66% leads to an asymme- 
try of the ti background of the same size as that of the 
signal. For the special case chosen with high px lep- 
tons in the final state, i.e. the high-py scenario with 
A231 = 0.3, a simulated asymmetry of 0.89% would 
lead to the same effect. Therefore, a higher pT cut is 



less sensitive to systematic errors, due to the high px 
cut effecting the tt background. 



V. CONCLUSION 

For the special case of dominant B3 couplings A'^gj, 
resonant charged slepton production is not possible at 
hadron colliders, as there are no incoming top quarks 
in the proton. We must then consider the associated 
production with a top quark. We have analyzed this 
difficult signature in detail. As the decisive observ- 
able for X[^i and A231 we propose the lepton charge 
asymmetry. For the supersymmetric spectrum SPSla' 
we have found a significant sensitivity range at the 
LHC, summarized in Fig. [TOl For the heavier spec- 
trum SPSlb the LHC is significantly less sensitive in 
the coupling, as can be seen in Fig. [11] We have 
constructed a heavy spectrum with a larger slepton- 
neutralino mass difference in order to explore the max- 
imum sensitivity at the LHC. This is shown in Fig. \T% 
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